A macro-realism inequality for opto-electro-mechanical systems 
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We show how to apply the Leggett-Garg inequality to opto-electro-mechanical systems near their 
quantum ground state. We find that by using a dichotomic quantum non-demolition measurement 
(via, e.g., an additional circuit-QED measurement device) either on the cavity or on the nano- 
mechanical system itself, the Leggett-Garg inequality is violated. We argue that only measurements 
on the mechanical system itself give a truly unambigous violation of the Leggett-Garg inequality 
for the mechanical system. In this case, a violation of the Leggett-Garg inequality indicates physics 
beyond that of "macroscopic realism" is occurring in the mechanical system. Finally, we discuss the 
difficulties in using unbound non-dichotomic observables with the Leggett-Garg inequality. 
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They [l8j showed side-band cooling from photo-pressure, 
and evidence of normal- mode splitting, i.e., strong cou- 
pling between optical and mechanical modes. Recent re- 
sults |15l Il6| using opto-electro-mechanical systems (i.e., 
a microwave transmission line in place of the optical cav- 
ity) have shown ultra-strong coupling and ground-state 
cooling. However one cannot easily distinguish the re- 
sultant effective low-temperature state of two coupled 
quantum oscillators from two coupled classical oscillators 
[l9l l2(ij . It is well known from quantum optics that the 
linear response spectral pro perties one observes are simi- 
lar for both theories [l8l l2ll. l22j , though spectral proper- 
ties can stron gly infer cooling to the mechanical ground 
state pE [la Il9l. [23j . In addition, the observation of 
asymmetry between spectral peaks due to absorption and 
emission of quanta is purely a quantum effect [22j, and 
has been recently observed in experiment 24 1 . 

In this work we propose a method of further dis- 
tinguishing quantum and classical oscillators by apply- 
ing the Leggett-Garg inequality. By using dichotomic 
quantum non-demolition measurements (QND) [IE [26[ 
we show that a theoretical model of a realistic opto- 
mechanical system implies a violation of the Leggett- 
Garg inequality due to the coherent interaction of the 
cavity with the mechanical oscillator. We show how ei- 
ther measurements on the cavity, or on the mechanical 
system directly, produce violations of the inequality. We 
argue that the latter are stronger proofs of quantum be- 
haviour in the mechanical system, as the former can also 
occur due to the quantum nature of the cavity alone. 

Since the dichotomic QND measurements we use here 
require strong coupling to a qubit our results are most 
directly applicable to opto-electro-mechanical systems 
which employ microwave transmission lines (as the "opto- 
electro-" cavity that cools the mechanical system to 
its quantum ground state). As far as we are aware, 
a dichotomic number state measurement has not been 
achieved in optical cavities. We believe the results we 
show here align well with Leggett and Garg's original 
goal of testing for non-realism in the macroscopic world, 
since the ground state (or a single Fock state) of a me- 
chanical oscillator represents a quantum state in a solid 



The Leggett-Garg (LG) inequality [l-|8j is one of a 
large class of inequalities used to delineate different phys- 
ical theories. It is constructed to test for "macroscopic 
realism" , the class of physical theories that imply that 
before we measure a property of a system, that property 
has a well defined value (which is not the case in quan- 
tum mechanics). Bell's inequality [1[ also tested for this 
property, but not without also testing for non-locality. 
The LG inequality [I|, attempts to test only for real- 
ism, but to do so requires the assumption of non-invasive 
measurement, and macroscopically distinct states. Hence 
the moniker of "macroscopic realism" . 

In the original LG proposal [I[ , they imagined measur- 
ing the two different and distinct "macroscopic states" 
of a superconducting flux qubit (where, mathematically, 
one can describe these states as a quantum two-level sys- 
tem). However, physically these two states are by most 
definitions "macroscopic" : They involve millions of par- 
ticles. Superconducting qubits have been used to show 
violations of Bell's inequality [I(|, the Leggett-Garg in- 
equality [5] , and have been proposed as a way to test the 
Kochen-Specker theorem [111 ]. 

An alternative candidate to test for quantum behavior 
in the macroscopic limit is in the ground state of a nano- 
mechanical oscillator [l2|, [l| . Strong evidence has been 
reported of success in this goal by coupling a nanome- 
chanical resonator to a qubit 14] . Recent work suggests 
that the ground state has also been reached in an opto- 
mechanical device [IE [I|| . An opto- mechanical system 
is essentially an optical (or microwave) cavity coupled to 
a mechanical resonator to cool and measure the mechan- 
ical system[T3, H3- A generic physical model for this 
opto-mechanical system is of a spring that supports one 
of the mirrors of an optical cavity, and thus the mechan- 
ical motion of the spring is coupled to the frequency of 
the optical mode. However, the physical realization of 
opto-mechanical devices can vary greatly, from a mirror 
suspended on a cantilever 18], to a mechanical mem- 
brane capacitively coupled to a microwave transmission 
line [IE El- 

Reference [H[ is an interesting example of the 
optical-cavity realization of an opto-mechanical system. 
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composed of millions of atoms. 

We begin this article by outlining the original Leggett- 
Garg inequality, and discuss why dichotomic QND mea- 
surements are necessary. We then show our main re- 
sult: that the introduction of a single photon into the 
microwave cavity, and application of dichotomic QND 
measurements, leads to a violation of the LG inequality. 
Afterwards we present the technical details of the model 
we use to describe the opto-electro-mechanical system, 
and discuss the practical issues of state preparation and 
measurement. We finish with a discussion of the difficul- 
ties of using non-dichotomic unbound measurements, and 
give a conjecture on a possible bound for the inequality 
in such a case. 



I. THE LEGGETT-GARG INEQUALITY 

The Leggett-Garg inequality [f], [|| is defined as fol- 
lows: Given an observable Q{t), which is bound above 
and below 

0,3 by \Q(t)\ < 1, the assumption of: (Al) 
macroscopic realism and (A2) non-invasive measurement 
implies, 



L(h,t 2 ) = (Q(h)Q(0)) + (Q{h- 
- (Q(h+t 2 )Q(0)) < l, 



-* 2 )Q(ti)> 



(i) 



If Q is in the steady state at the initial time of measure- 
ment, and we set t\ = t 2 , then this becomes 



2(Q(r)Q(0)> - <Q(2t)Q(0)> < 1. 



(2) 



To adapt this to work on measurements on bosonic 
(harmonic) systems one must proceed with extreme cau- 
tion. This is because (a) it is difficult to define a bound 
on measurements on harmonic systems, (b) many mea- 
surements (particularly in the optical regime) are inva- 
sive (e.g., single-photon counting) and (c) the dynamics 
of classical and quantum harmonic systems are identi- 
cal (apart from quantum fluctuations) without additional 
sources of non-linearity. 

Fortunately, the growing field of op to- 
electromechanical systems and circuit QED |27H3lj 
allows us to overcome many of these obstacles. We 
can adapt the scheme realised by Johnson et al [l^, [2(| 
to overcome obstacle (a). In their scheme, one uses 
an additional qubit/measurement-cavity system to 
dispersively measure whether the optomechanical-cavity 
contains "one photon or not" (this is a dichotomic QND 
measurement). We will also show how, in principle, it 
might be possible to use this to measure the mechanical 
system directly, and say if it contains "one phonon or 
not". 

The former (a dichotomic QND measurement on the 
cavity) is possible due to the strong coupling between 
qubit and cavity that has been achieved in circuit-QED 
systems. The latter (a dichotomic QND measurement on 
the mechancical system) may be possible given the recent 
strong coupling shown between a mechanical system and 



a superconducting qubit [Tj|- A possible realization of 
the QND measurement on the cavity in an opto-electro- 
mechanical system is shown in Fig. 1. 

This scheme also allows us to overcome obstacles (b) 
and (c), as it realizes a non-demolition, and classically 
non-invasive, projective measurement of the photons in 
the cavity (or phonons in the mechanical oscillator). In 
the final section we will return to these issues, and con- 
jecture about a new bound for the inequality if one's 
observables are non-dichotomic and unbound. 
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FIG. 1: (Color online) Schematic diagram of a possible 
opto-electromechanical system with dichotomic QND mea- 
surement of the cavity. Right side of the figure: the driven 
microwave cavity is coupled to a mechanical oscillator, e.g . 
the fundamental mode of a thin film drum as in [l5|, Il6{ . 
The left side of the figure shows a simplified schematic of 
the qubit/measurement-cavity system used for the dichotomic 
QND readout (as in Ref. 25]). A similar measurement could 
be performed directly on the phonon states in the mechan- 
ical mode, given sufficiently strong coupling strengths be- 
tween qubit and mechanical system (such a configuration is 
not shown in the figure). 



II. VIOLATION FOR OPTOMECHANICAL 
SYSTEMS USING DICHOTOMIC QND 
MEASUREMENTS 



We define the LG inequality in terms of dichotomic 
quantum non-demolition measurements either on the 
single-Fock state occupation of the cavity mode 



2|l)c(l| 



1, 



(3) 



where c refers to the cavity mode, or on the single-phonon 
state occupation of the mechanical mode, 



211) 



(II 



1, 



(4) 



where m refers to the mechanical mode. As men- 
tioned above these measurements requires an additional 
qubit/measurement-cavity (25j . which we outline in sec- 
tion IIV1 and is shown schematically in Fig. 1, for the 
example of measuring the cavity mode. For our pur- 
poses, this measurement returns +1 if there is a quanta 
in the appropriate mode (m or c), and —1 if not. To 
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show a violation of Eq. ([T]) one can prepare the opto- 
electromechanical system near its ground state following 
the side-band cooling procedure described in the next 
section. The ground state cooling of the mechanical sys- 
tem requires that we strongly drive the microwave cavity, 
resulting in a non-zero steady-state coherent occupation 
in the cavity in a rotating frame. Fortunately this can, 
in principle, be eliminated from the QND measurement 
(see section HVj) . 

We then adiabatically introduce an additional photon 
into the cavity [32j , in addition to this coherent state [H, 
HI], which ideally prepares the system in the state 

p(t = 0) = |l) c (l| c ®|0) m (0| m , (5) 

where again c refers to the cavity and m refers to the 
mechanical mode (note that the state of the cavity is 
in a displaced basis in a rotating frame because of the 
driving of the cavity). The strong coupling between the 
mechanical system and the optical mode causes this sin- 
gle excitation to be coherently exchanged (akin to a Rabi 
oscillation). One then measures the operator Q c (or Q m ) 
using the readout qubit-measurement-cavity and a pro- 
grammable C-NOT scheme [2|| (see section HV|) . If the 
measurement timescale (which includes rapidly resetting 
the qubit to its ground state) is short enough one can 
construct the two-time correlation functions in Eq. (fl}. 

In Fig. 2 we explicitly show how the results from our 
model, outlined in the next section, which suggest a vio- 
lation of Eq. Q, is in principle observable with existing 
experiments (14H16| . It is interesting to note that the 
largest violations occur for small times, which implies 
we require the readout and reset of the qubit to be fast. 
Typical readout times in Ref. [25[ are of the order of 
550 ns, introducing an intrinsic minimum delay into the 
correlation functions. The typical time scale of the coher- 
ent dynamics in Ref. [H, [Ty] is related to the coupling, 
which is of the order of 10 5 -10 6 Hz, implying that the 
short time-scale coherent dynamics should be observable 
with such a measurement, though for measurements of 
Q c this may be altered by the need to measure in a ro- 
tating frame (see later). We also find that if the driving, 
and hence the coupling between cavity and mechanical 
system, is strong enough then the non-energy conserving 
terms in the interaction modulate the dynamics quite 
strongly. However at this point one also expects other 
non-linear affects to arise. 



A. Ambiguities in cavity measurements 

As we will discuss shortly, the measurement of Q c is a 
direct adaptation of an existing experiment (albeit with 
additional steps to make sure that we are measuring in 
the correct frame). However, our main goal is to ver- 
ify the quantum dynamics of the mechanical mode. It 
just so happens that in this case it is the quantum coher- 
ent interaction between the cavity and mechanical modes 
which drives the violation we observe in the observables 



of the cavity system. However, in principle a violation 
could also be observed due to the quantum nature of the 
cavity mode alone. 

Thus with measurements on the cavity mode alone it 
is impossible for us to state that a violation the Leggett- 
Garg inequality (e.g., with measurements Q c ) gives un- 
ambiguous proof of macroscopic quantum phenomena in 
the mechanical mode. Ideally, one requires dichotomic 
QND measurements on the mechanical mode directly (as 
defined by Q m ) to state that a violation of the Leggett- 
Garg inequality is unambiguous proof of quantum me- 
chanics in the nano- mechanical system. We will outline 
a possible scheme to achieve this later. 

III. OPTOMECHANICAL SYSTEMS 

We now explicitly describe the optomechanical system, 
and the model we use to calculate the results shown in 
Fig. 1. This model is well known and studied in other 
works [H, [34| , but we provide details here for clarity. 

We start with the Hamiltonian describing the coupling 
between the cavity and the mechanical oscillator |35| . 

H (1) = Aa^ a + ui m tf b + gu> m (b + tf)a^a 

+ fi(a + a t ) (6) 

The driving field is such that the cavity and mechanical 
mode are now near resonance (A = u) m ). 

One of the approaches taken before (e.g., dHHH) is 
to insert displacements a and /3 for both modes, take the 
limit g,ga <C A,w m , and treat the cavity mode as an 
effective environment which cools the mechanical mode. 
The condition for cooling, in our notation, is then A > 0, 
and a sufficiently high-quality cavity k < uj m . Reference 
[H ] has shown that by observing the homodyne transmis- 
sion spectra from the (optical) cavity, we can see a clear 
signature of "normal mode splitting" . This is because the 
displacements a and /?, and the coupling, are given by a 
simple model of two coupled oscillators, which has cou- 
pled normal modes. Here we explicitly model both cavity 
and mechanical system in the strong-coupling regime ob- 
served in Ref. [151 ] using a master equation approach. 

A. Resolved side-band cooling 

One can find the linearized version of the original 
Hamiltonian by displacing both mechanical and cavity 
modes, so that a — > c + a, b — s- d + (3. Inserting these 
displacements in the Hamiltonian, and eliminating linear 
terms, gives two coupled equations for the displacements, 

aA + 2w m ga{P + p*) + fl - in = (7) 

oJm/3 + 0J m ga 2 = (8) 

The k term arises because of the linear dissipation 
terms (cavity losses) we will introduce shortly. Careful 
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inspection of the possible solutions of these cubic equa- 
tions shows that, 



is to prepare the system in the appropriate initial state; 
e.g. a single photon in the cavity (in the displaced basis) 



P = —gaa*, 



and in the limit of small g, 



-n 



A - in/2' 



(9) 



(10) 



If we do not make a small g assumption, these displace- 
ments are real, up to some critical driving Q of order 
Um/g, which corresponds to the breakdown point of the 
small displacement assumption made to derive the origi- 
nal Hamiltonian (see Ref. [35[). At this point additional 
non-linearities in the interaction could play a role, but 
we do not consider those here. 

The Hamiltonian, with the linear terms eliminated, be- 
comes 



H {2) = (A + 2guj m /3)c' !< c + uj m dU 



(11) 



+ guj m (d + d')(a*c + ac^) + guj m (d + d')c'c. 

We then add standard Lindblad cavity and mechani- 
cal losses to this model, and solve the resulting Master 
equation, 

p = -i[H™,p] 

+ — {— c'cp — pc^c + 2cpc* 

+ \(ot*c — ac^p + p(ac ! — a*c)l } 

+ 2^ + r > [- dfd P ~ P dtd + 2d P d ^] 
+ [-dSp - pdd^ + 2<tfpd] 

+ ~ [((3*d-0)p + p(f3j ~/3*d)] 

(12) 

Where N is the initial thermal occupation of the me- 
chanical mode. Dissipation terms linear in c and d (and 
displacements a, /3) arise because of the shifted coordi- 
nate frame. As mentioned earlier, the linear terms for 
the cavity can be easily eliminated by including them in 
the displacement a. The linear terms for the mechanical 
mode dissipation are small in the limit of a high qual- 
ity factor resonator, so we neglect them here (though we 
have numerically checked that their influence is small). 

Under the conditions k < u> m , A > and sufficiently 
large driving strength Q, one can achieve the well-known 
resolved side-band-limit cooling; one can start from a 
thermal state of the resonator at a given temperature, 
and reach a steady-state, where the thermal phonon oc- 
cupation of the mechanical system approaches zero. See 
Refs. [Ill, [HJ for further details and discussion of the 
cooling process. 

Using this model we can easily construct the various 
correlation functions needed for Eq. (TTJ). Our technique 



p(0) = |l C) m )(l c ,0 m |, 



(13) 



then the appropriate correlation functions are calculated 
via the time evolution 

(Qi(2T)Q i {r))= r Tr[Q i exp[LT}Q i ex V [LT]p(p)], (14) 

or via the quantum regression theorem. 

Since we operate always in the basis of the displaced 
modes, we are always close to the steady state. Thus 
imposing directly as initial conditions a single Fock state 
is a sufficiently good approximation to the true process 
of preparing the opto-mechanical system in its steady 
state, and then, e.g., introducing the single-photon state 
using the measurement qubit [25|. In principle one can 
explicitly model this state-preparation stage [33J , but for 
simplicity we omit it here. 



IV. QND READOUT 

As discussed earlier, both of these measurements, Q c 
and Qrm are challenging, but may be feasible in the 
future by combining existing circuit-QED devices (for 
QND readout [25[) with an opto-electro- mechanical sys- 
tem [l5|, [l6| . The additional circuit-QED system (qubit 
and microwave cavity) allows both the deterministic 
preparation of the cavity in a single Fock state [26|, |36| - 
|38| |. and the dispersive QND readout of its population 
dynamics [25[ . Thus in reality our proposed opto-electro- 
mechanical system is a circuit-QED-mechanical system, 
where the additional qubit-cavity part is used for state 
preparation and readout. First of all we will describe 
the details (Fig. 1) of how to realize the measurement 
Q c (Ref. [25(). As mentioned earlier, this is perhaps the 
most feasible with current technology, though does not 
give us an unambigous violation of the LG inequality for 
the mechanical mode. Then we will discuss possible ways 
that the measurement Q m might be realised by coupling 
the qubit directly to the mechanical system, and not the 
cavity, which give us a more ideal and unambigous vio- 
lation of the LG inequality. 



A. Cavity measurement Q c 

The dichotomic QND measurement realized by John- 
son et al [25[ is ideal for our purposes of realizing Q c , but 
the scheme as it is described there [25[ measures "if there 
is one cavity or not" in the lab basis. The cavity in the 
opto-mechanical Hamiltonian Eq. © we used earlier is 
in a displaced rotating frame (because of the microwave 
driving needed for sideband cooling), and thus it is in 
this basis that we must measure the cavity to realize Q c 
as we have described it. For example, in the stationary 
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frame of the qubit (but rotating frame of the cavity), the 
interaction between qubit and cavity is described by, 



H q b-c — 2°" z + Aata + A (er_| 



(15) 



where LUd is the driving frequency, and was chosen to 
bring the cavity and mechanical system on resonance in 
Eq. ©, so that 



(16) 



In addition, we also displace the cavity co-ordinates by 
a, so that the interaction between the qubit and cavity 
co-ordinates that we actually want to measure is 



H, 



qb-c = -Z a z + Ac'c+ A (<T + Ce 1- rj r'< 



-iu d t 



+ A (ae~ iUdt (T + + 



a e <j-\ 



(17) 



The additional displacement term represents the large 
number of photons that are in the cavity due to the driv- 
ing. Ideally their influence on the qubit can be eliminated 
by ap ply ing an additional microwave drive to the qubit 
itself [28[ (still in the lab frame), out of phase with the 
term above, e.g., 



H, 



adjust 



= -A 



-iuidt. 



(18) 



This is feasible if the magnitude, Act, is not too large [28j j . 
but may become unfeasible if an extremely large driving 
of the cavity is needed for cooling. 

Assuming this term has been applied, and the effect of 
the large cavity population eliminated, we can move the 
qubit into the same frame as the cavity with the unitary 
transformation U q = exp (iujd<J z t/2). This leaves us with 
a normal Jaynes-Cummings Hamiltonian between qubit 
and cavity, with a shifted qubit energy 



A' 



(19) 



shift. In addition, the higher-order terms in X/5 (repre- 
senting back-action of the qubit on the cavity) should be 
much smaller than the cavity-mechanical mode interac- 
tion (i.e., A 3 /(5 2 <C gaoj m ). 

In Ref. (25| , in order to have sufficiently high resolution 
measurement of the effect of the photons on the energy 
levels of the qubit, they needed 5/ A < 10. That is, they 
need sufficiently large S to reach the dispersive limit, but 
sufficiently strong A to obtain well resolved energy shifts 
for different photon occupations. Here, S == e — uj c , thus 
reaching the same regime as <5 / A < 10 seems feasible. 

The dichotomic property of the measurement is 
achieved because of the strong dependence of the qubit 
response on the number of photons in the cavity. In 
Ref. [25j . for the measurement step, they apply a ir con- 
trol pulse to the qubit at the frequency corresponding to 
its energy when just one photon is in the cavity. Thus 
the qubit is rotated if and only if there is one photon 
present, and nothing happens otherwise. This is an ef- 
fective CNOT gate on the qubit and the cavity. Here, the 
effective CNOT gate must also be in the rotating frame 
of the qubit. For the final measurement step, one mea- 
sures the state of the qubit via pulsed spectroscopy of 
the cavity. 

Overall, this dispersive Hamiltonian, combined with 
the controlled- 7r-rotation of the qubit, and readout of the 
qubit using the additional measurement cavity (which 
we have not explicitly described), ideally gives us a 
way to realize the dichotomic QND measurement Q c — 
2|l)c(l| e — 1. As we discussed earlier, the time needed in 
Ref. [251 ] to realize this measurement may be short enough 
to observe correlation functions on the time scale we re- 
quire. However, in general there will be losses involved 
in the measurement process (e.g., due to dissipation of 
the qubit state) which will degrade the measurement re- 
sult [25|] . In addition, the need to perform the effective 
CNOT gate in the rotating frame may slow down the 
measurement step (28|. 



In this new picture, the QND measurement scheme 
proposed and analyzed elsewhere [HI, [H, [3!| applies for 
a large bias 



5 = A' - A 



(20) 



This is clearly shown by applying the unitary transfor- 
mation UdHU^, Ud = exp [-|(cct+ — ct<r_)], which leads 
to the well-known dispersive coupling Hamiltonian (see, 
e.g., Refs. [HHI), 



H, 



qb-D 



T 



A 2 

5 



CC+ 2 



a z + Ac f c. (21) 



This transformation can induce interactions between the 
measurement qubit and the mechanical mode, but these 
terms can also be treated with a dispersive transforma- 
tion, and give a shift of the qubit frequency of order 
(Xgauj m ) 2 / S 3 , and are thus much weaker than the A 2 /<5 



B. Mechanical measurement Q m 

As we have reiterated several times, a measurement 
of Q c is not sufficient to unambiguously show quantum 
dynamics in the mechanical mode. We ideally need to 
perform the dichotomic QND measurement on the me- 
chanical system itself. As far as we are aware no sim- 
ilar measurement has yet been achieved, though efforts 
on membrane-in-the-middle devices [22} are promising. 
Staying within the regime of the nano-mechanical sys- 
tems we have discussed so far [I3 - [T^I. one can imagine 
adapting the scheme of Johnson et al 1251 to directly mea- 
sure the mechanical mode. In Ref. [Ij] O'Connell et al 
observed a strong interaction between a high-frequency 
mechanical mode and a superconducting qubit. There 
the mechanical mode frequency was w m /27r = 6 GHz, 
so cryogenic freezing was sufficient to reach the quantum 
ground state, and they observed qubit-oscillator coupling 
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strengths of X/h = 110 MHz. This is favorable for using 
the qubit as a dispersive measurement of the mechan- 
ical mode. However a straight adaptation of Ref. [25[ 
to the system in Ref. [3] would have to compensate for 
the extremely short quality factor of the mechanical res- 
onator (the mechanical dephasing time is estimated to 
be T w 20 n s). This is well short of the measurement 
time in Ref. [25| [. and thus resolving the short time cor- 
relations needed to see a violation of the inequality may 
prove difficult without improvements in the readout and 
reset times of the qubit, or employing a higher quality 
factor /lower frequency mechanical resonator. 

Furthermore, one can imagine a similar scenario using 
the opto-mechanical side-band cooling systems we have 
described here, where the low frequency (and high qual- 
ity factor) mechanical system is cooled by the cavity, and 
then the mechanical part is measured in the same man- 
ner as above (by an additional superconducting qubit, 
with compensation for the coherent occupation 0). This 
is quite a speculative scenario, as it is not clear if a suf- 
ficiently large coupling between the qubit and the type 
of mechanical oscillator used in opto-electro-mechanical 
devices [H, [l6[ can be engineered, and if the overly large 
energy mismatch between the superconducting qubit and 
mechanical resonator overcome. However if realised it 
would be ideal for truly showing macroscopic quantum 
phenomena in the same spirit as Leggett and Garg's pro- 
posal [lj. 

C. Single photon measurements 

For the optical-cavity realization of an optomechancial 
system (as in Ref. [Hj]), the Fock state preparation (e.g., 
by u sing an additional one-way cavity as discussed in 
Ref. [33]), and QND measurements [40] , are feasible but 
the dichotomic measurements we require are much more 
difficult to realise than in the microwave cavity case. 

We also point out that the correlators in Eq. (fTJ) are 
not normal ordered, and thus do not represent the mea- 
surements obtained from single-photon counting (which 
are typical in optical cavity systems). As we discussed 
in earlier work [6|, @] , photon absorbtion measurements 
are fundamentally invasive, and typically represent an 
obstacle for un-conditionally verifying quantum behavior 
via the Leggett-Garg inequality. This is particularly true 
with a fragile single-quantum Fock state, hence the need 
for QND measurements. 



V. NON-DICHOTOMIC AND UNBOUND 
OBSERVABLES 

What happens if we attempt to construct the Leggett- 
Garg inequality from non-dichotomic and unbound ob- 
servables? Recent work on Bell's inequality with un- 
bound measurements fill . [42j suggest that one has to 
move to fourth-order correlation functions to distinguish 



quantum and classical correlations, which may also apply 
to the Leggett-Garg inequality. 

First of all, let us consider a general picture where 
we measure an unbound operator (Q) G {— oo,oo}. Fol- 
lowing the same reasoning as used in the Leggett-Garg 
inequality one can derive a bound (and assuming we con- 
struct our expectation values by counting how often a 
particular measurement result arises), 

L Q < (Maxi[Q(t)] 2 }. (22) 

Such a bound may occur due to some intrinsic conserva- 
tion rule in the system (e.g., if the number state or energy 
is conserved). However, this bound is both difficult to 
calculate (and measure) and is extremely loose, since in 
general extremum values may be observed but contribute 
little to the expectation values. Since the maximization 
is a convex function, we know that 

Max t [(Q(i) 2 )] < (Max t [Q(t)] 2 ), (23) 

(this is the Jensen inequality), but finding further con- 
straints on these functions is challenging beyond trivial 
cases. We conjecture that there might be a tighter bound 
for the inequality given by Max t [(Q(i) 2 )]. However we 
have been unable to find a rigorous proof, and it may be 
that a simple counter-example exists to show that this 
conjecture does not hold. 

There is a further caveat on such an approach. Note 
that taking the stationary limit and setting (ti — t2 — 
t) simplifies the inequality with our conjectured bound 
Max t [(Q(i) 2 }] to, ig' = 2(Q(r)Q(0)) - (Q(2r)Q(0)> < 
(Q(0) 2 }. In a classical situation, the correlation functions 
one can observe are harmonic functions (even in the sta- 
tionary state). For example, one can easily solve the 
equation of motion for a single oscillator in contact with 
a thermal bath, and find that the spectral density of the 
displacement [12| is a Lorentzian with amplitude re- 
lated to the bath temperature T. The Wiener-Khinchin 
theorem tells us that the spectral density is the Fourier 
transform of the auto-correlation function in the steady 
state, implying sinusoidal correlation functions which are 
only dependent on one time variable (the time between 
measurements). These obviously cause a violation of the 
bound (Q(0) 2 } in the steady state. One can argue that 
this is not per se a failure of our conjectured bound, but 
is because the only output one observes from the sys- 
tem is noise-driven; the thermal background has a white 
noise spectrum, which can thus excite the system around 
its resonant frequency. Thus, in the language of Leggett- 
Garg, these classical correlation functions are essentially 
invasive since one only observes a signal when the system 
fluctuates (e.g., due to thermal fluctuations). Therefore, 
the observation of thermal-noise-induced fluctuations is 
equivalent to a perturbation of the system by the mea- 
surement. 

The original Leggett-Garg inequality avoids this prob- 
lem by demanding that the system must be in one of 
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two macroscopically distinct states, and that the sys- 
tem is almost always in one of the two states [l|, 
In a harmonic system this assumption, of macroscopi- 
cally distinct states, breaks down spectacularly. One can 
overcome this problem by avoiding the steady state, or 
introducing a third measurement (if Q G {0, oo}, or a 
third and fourth measurement if Q G {—00,00}) into all 
correlation functions in the inequality at time t = 0, and 
scaling the bound appropriately. Then the violation is de- 
pendent on the effect of the second measurement (which 
we assume again to be non- invasive) . The quasi- invasive 
(fluctuation) nature of the first measurement becomes 
irrelevant. Introducing extra measurements into the in- 
equality is akin to the fourth-order Bell inequality derived 
by Bednorz et al (4lL |42j ■ However more work remains 
to be done to derive rigorous proofs for our conjecture. 

VI. CONCLUSIONS 

In summary, we discussed how to use the Leggett-Garg 
inequality to distinguish quantum and classical dynamics 
in opto-electro-mechanical systems. We illustrated that 
dichotomic QND measurements of either the cavity or 



mechanical system leads to a violation of this inequality. 
We discussed possible methods to realize such measure- 
ments, and argued that only measurements directly on 
the mechanical system itself will give unambigous proof 
of macroscopic quantum dynamics in the mechanical sys- 
tem (or, in the language of Leggett and Garg, proof of a 
violation of macroscopic realism). 
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FIG. 2: (Color online). Example of violation of Eq. ([T]) 
(with ii = £2 = r) using (a) measurements on the cavity 
Q c and (b) measurements on the mechanical mode Q m , as 
a function of the dimensionless time: t\G\/2h. Both figures 
are produced by our model of an opto-electromechanical sys- 
tem, where G = gui m ct is the effective coupling between the 
mechanical and cavity modes, and a ~ — S7/A is the displace- 
ment of the cavity mode produced by the driving SI. For 
strong effective coupling, the non-energy conserving terms in 
our model begin to strongly modulate the shape of the cor- 
relation functions in Eq. ((TJ). We have chosen ratios for the 
parameters that correspond ap pro ximately to those seen in 
the strong-driving limit of Ref. [15|, [l6| , where u m ~ 2ir x f 
MHz, though to aid computation we made g large, and Q 
small, in comparison to their data. Remarkably, the high 
quality factor of both the microwave and mechanical cavities 
used in [HI . [Tr| means that the violation of the inequality re- 
mains visible for relatively long time scales, though this can 
depend on the initial temperature of the mechanical mode. 



